PARTON DISTRIBUTIONS IN NUCLEI AT SMALL X 
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We study the nuclear shadowing effect in the context of Glauber-Gribov multiple- 
scattering model and perturbative QCD. We find that at small x, the Q 2 evolution 
of the shadowing is much slower than the DGLAP evolution, due to the multiple 
scatterings at small x. We show that the gluon shadowing at small x and for 
Q 2 > 3GeV 2 is perturbative in nature and does not depend on the initial, non- 
perturbative condition. We evaluate the impact parameter dependence of the gluon 
distribution and show that it is a non-linear effect in the nuclear thickness function. 

The semihard processes, such as minijets, dileptons, open charm and di- 
rect photon production in heavy-ion collisions at y/s > 200 GeV, can be reli- 
ably calculated in the framework of perturbative QCD assuming the validity 
of the factorization theorem in perturbation theory, and with the knowledge 
of parton distributions in nuclei at small x. 

Quark distribution in nuclei has been measured and the attenuation has 
been firmly established experimentally at CERN El and Fermilab □ in the region 
of small x in deeply-inelastic lepton scatterings (DIS) off nuclei. The data, 
taken over a wide kinematic range, 10~ 5 < x < 0.1 and 0.05 GeV 2 < Q 2 < 100 
GeV 2 , show a systematic reduction of nuclear structure function F% (%, Q 2 )/A 
with respect to the free nucleon structure function F% (x, Q 2 ). However, there 
is no experimental measurement of the gluon density in nuclei. 

At low Q 2 , in DIS, the interaction of the virtual photon with the nucleons 
in the rest frame of the target is most naturally described by a vector-meson- 
dominance (VMD) model 0. At Q 2 > 1 ~ 2 GeV 2 , the virtual photon can 
penetrate the nucleon and probe the partonic degrees of freedom where a 
partonic interpretation based on perturbative QCD is most relevant in the 
infinite momentum frame. In the target rest frame, the virtual photon inter- 
acts with nucleons via its quark-antiquark pair (qq) color-singlet fluctuation 
□. If the coherence length of the virtual photon is larger than the distance 
between nucleons in a nucleus, l c > Rnn, the QQ configuration interacts co- 
herently with fraction of the nucleons, while for l c > Ra (i.e. x < 10~ 2 ), 
it interacts coherently with all the nucleons, with a cross section given by 
the color transparency mechanism for a point-like color-singlet configuration, 
PqqN = ^-r 2 a s x' guLA^x' , 1/r 2 ) i. This results in a reduction of the total 
cross section and consequently attenuation of the parton distributions in nu- 
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In the Glauber-Gribov multiple-scattering theory U nuclear collision is a 
succession of collisions of the probe with individual nucleons within nucleus. 
A partonic system (h), being a qq or a gg fluctuation, can scatter coherently 
from several or all nucleons during its passage through the target nucleus. The 
interference between the multiple scattering amplitudes causes a reduction of 
the hA cross section compared to the naive scaling result of A times the 
respective hN cross section, the origin of the nuclear shadowing. The total 
hA cross section is given by 



Jd 2 h2 [l - e-'w^Cbyaj =2irR 2 A [ 1 + ln Kh + E 1 ( Kh )} , (1) 



where Kh — AahN / {^R\) is an impact parameter averaged effective number 
of scatterings. For small value of Kh, o~hA — * 2irR 2 A Kh = Ao~hN, the total 
hA cross section is proportional to A. In the limit Kh — > oo, the destructive 
interference between multiple scattering amplitudes reduces the cross section, 
o~ hA — * 2nR A (-f + hi Kh)- Namely, the effective number of scatterings is large 
and the total cross section approaches the geometric limit 2nR A , a surface 
term which varies roughly as A 2 / 3 . 

In the Glauber-Gribov eikonal approximation, 



o 



o-( 7 *A)= d 2 b dz d 2 r\ij{z,r)\ 2 2(l-e-^ N ^^ T ^ b ^ 2 ), (2) 



where |^i(j»,r)| is the photon wave function and z is the fraction of the 
energy carried by the quark (antiquark). The nuclear cross section is therefore 
reduced when compared to the simple addition of free nucleon cross sections. 

At small x, the structure function of a nucleus, F^{x, Q 2 ) can be obtained 
from a(-y*A). Substituting integration over (z, r t ) to (x',Q' 2 ) in (2), one 
obtains for (sea)quarks □ 

x f A (x,Q 2 )=xf A (x,Q 2 ) + ^4^ 4/ dQ«[ 7 + ln(«,)+J5i(« g )](3) 

07T J x X Jq2 

where K q (x,Q 2 ) = 3 ^fg 2 a s {Q 2 )xg^ LA (x, Q 2 ). We find our result for 

F^(x,Q 2 ) to be in excellent agreement with thej-ecent HERA data in the 
kinematic region of x and Q 2 where DLA is valid B. 

We parametrize the initial shadowing ratio Ro(x) = 
F^(x,QI)/AF^(x,Qq) at Ql = 0.4 GcV 2 using experimental data and 
we calculate the nuclear structure function at the measured (Q 2 ) values at 
different x values. The results for 40 Ca and 208 Pb are shown in Fig. 1. 
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We note that at large Q 2 the nuclear shadowing effect is reduced but does 
not diminish, i.e. it is not a higher twist effect. This can be understood as 
the interplay between the perturbative and the non-perturbative shadowing 
mechanisms, which is also evident in the ^-dependence of R q El. The apparent 
flatness of the shadowing ratio at low Q 2 in the small- a; region is altered by 
the perturbative evolution. This is due to the singular behavior of xg® hA 
as x — ► at large Q 2 leading to the strong x-dependence of the effective 
number of scatterings, K q (x). Furthermore, at small x the Q 2 -evolution of 
the F2/AF2 is much slower than the DGLAP evolution due to multiple 
scatterings. We illustrate this effect in Fig. 2. When the effective number 
of scatterings, n q (x), is small, the Eq. (3) becomes the DLA of the DGLAP 
evolution equation for quarks. 

Similarly, gluon distribution in a nucleus is given by 

2R 2 f 1 dx' f Q2 

xg A {x, Q 2 ) = xg A (x, Q 2 ) + — ^ / — / dQ' 2 [7 + ln( Kg ) + E\ {n g )\A) 

Jx % JqI 
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The essential difference between the quark and gluon cases is different splitting 
functions and the effective number of scatterings, K g = 9k 9 /4 due to different 
color representations that the quark and the gluon belong to. These two 
effects combined result in the 12 times faster increase of the gluon density 
with Q 2 than in the case of the sea quarks in the region of small x. The two 
important effects which make the gluon shadowing quite different from the 
quark shadowing are the stronger scaling violation in the semihard scale region 
and a larger perturbative shadowing effect. This can be seen by considering 
the shadowing ratio 



R g (x,Q 2 



xg A (x, Q 2 ) xg N {x, Ql)R° g (x) + Axg A {x; Q 2 



(5) 



Axg N {x,Q 2 ) xg N (x,Ql) + Axg N (x;Q 2 ,Ql) 

where R g (x) is the initial shadowing ratio at Qq and Axg(x; Q 2 , Qq) is 
the change of the gluon distribution as the scale changes from Qq to 
Q 2 . The strong scaling violation due to a larger n g at small x causes 
Axgiy(x; Q 2 7 Ql) ^> xg^(x, Qq) as Q 2 is greater than 1 ~ 2 GeV 2 . As seen in 
Fig. 3(a) for Q 2 > 3GeV 2 the dependence of R g (x, Q 2 ) on the initial condition 
R g {x) diminishes and the perturbative shadowing mechanism takes over. 

The x-dependence of the gluon shadowing can also be predicted as long 
as Q 2 > 3GeV 2 where the influence of the initial condition is minimal. The 
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shape of the distribution is quite robust in the small- a; region regardless of 
what initial conditions one may choose. Due to the perturbative nature of 
the shadowing, these distributions do not exhibit a saturation as x decreases. 
Furthermore since the shadowing is a non-linear effect in the effective num- 
ber of scatterings, the impact parameter dependent shadowing ratio cannot 
be factorized into a product of an average shadowing ratio and the nuclear 
thickness function. Our results are presented in Fig. 3(b). 

In summary, the nuclear shadowing phenomenon can be understood as 
a consequence of the parton coherent multiple scatterings. We find that the 
quark shadowing arises from an interplay between "soft" physics and the semi- 
hard QCD process. At small x, we find that the quark shadowing ratio evolves 
with Q 2 slower than the DGLAp evolution. We show that gluon shadowing 
is largely driven by a perturbative shadowing mechanism due to the strong 
scaling violation in the small- a; region. The gluon shadowing is thus a robust 
phenomenon at large Q 2 and can be unambiguously predicted by perturba- 
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tive QCD. The strong scaling violation of the nucleon structure function in 
the semihard momentum transfer region at small x can be reliably described 
by perturbative QCD and is a central key to the understanding of the scale de- 
pendence of the nuclear shadowing effect. The impact parameter dependence 
of gluon shadowing is a non-linear effect in the nuclear thickness function. It 
is important to correctly incorporate the impact parameter dependence of the 
nuclear structure function when one calculates the QCD processes in the cen- 
tral nuclear collisions. In the asymptotic limit of small x, gluon distribution 
in nuclei does not exhibit saturation, as in the case of the GLRlI equation, but 
rather tends to the limit xGa ~ R A Q 2 ln(l/x). Thus, the shadowing ratio, 
R g , does not saturate at small x, as long as the gluon density in a nucleon 
has singular x-dependence at small x. 
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